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Abstract 

Systems of hydrodynamic type equations derived from the Navier-Stokes 
equations and the boundary layer equations are considered. A transformation of 
the Crocco type reducing the equation order for the longitudinal velocity compo- 
nent is described. The issues of nonlinear stability of the obtained solutions are 
studied. It is found that a specific feature of many solutions of the Navier-Stokes 
equations is instability. The nonlinear instability of solutions is proved by a new 
exact method, which may be useful for the analysis of other nonlinear physical 
models and phenomena. 

Self-similar, invariant, partially invariant, generalized separable and some 
other exact solutions of the Navier-Stokes equations are considered, e.g., 
in [1-11]. 

Keywords: Navier-Stokes equations, boundary layer equations, Crocco trans- 
formation, exact solutions, nonlinear instability, nonlinear stability, Calogero equa- 
tion, exact methods. 

1 Systems of hydrodynamic type equations 

Systems of equations under consideration. The following systems of equation will 
be considered: 

g2p Q2p Qp 2 Q3p gp 

dG dG ,BF d^G 

— +F—-kG—=v— + s{t), (2) 
at oz az az-^ 

where F is one of the components of the fluid velocity, G is an auxiliary function, and 
V is the kinematic viscosity of the fluid. For m =1/2 and m = 1, a single equation 
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(1) and system of equations (1), (2) with m — 1 were treated in [5, 8, 10, 11], where 
exact solutions of the Navier-Stokes equations and the boundary layer equations were 
considered. Note that, for rri = 1, the functions p = p(t) and q = q{t) in equation 
(1) can be chosen arbitrarily [11]. Below, a new class of exact solutions of the Navier- 
Stokes equations will be described for system (1), (2) with m = k ~ 1/2 and q{t) = 0. 

Nonlinear equation (1) can be considered separately and equation (2) is linear with 
respect to the unknown function G. According to the general property of equation (1) 
[8, 1 1], if F{z, t) is some solution to this equation, then the function 

F^F{z + ij{t),t)-^P[it) (3) 

with an arbitrary function tp{t) will also be a solution to equation (1); in addition, the 
function F — —F{—z, t) will also be a solution to this equation. 

One class of exact solutions of the Navier-Stokes equations. Three-dimensional 
nonstationary motions of a viscous incompressible fluid are described by the following 
system of the Navier-Stokes and continuity equations: 

at ox oy oz 



p \ dx^ dy"^ dz 

^ + ^ + ^ ^0 
dx dy dz 

Here x, y, and z are the Cartesian coordinates; t is the time; Vi, V2, and V3 are the fluid 
velocity components; P is the pressure; p is the fluid density; and ViP = dP/dx, 
V2P — dP/dy, V3P — dP/dz. In equations (4), it is assumed that the mass forces 
are potential and included in the pressure. 

For a viscous incompressible fluid, equations (4) admit the following exact solu- 
tions: 

1 dF 1 dF 

^^-^-2^a^' ^^--2^a7' = 

- = Po{t) + ■^a{t){x^ + z/) _ s{t)x --F^ + j — dz, 

where pQ{t), a{t), and s{t) are arbitrary functions of time t and the functions F and G 
depend on z and t and satisfy the equations 

did-z^^^-A-d;) 

dG dG IdF d^G 

ar+^5^-2^a7-^a;^ + ^W. (6) 

System of equations (5), (6) is a special case of system (1), (2) with m — k — 1/2 and 
qit) = 0. 
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2 Reduction of order of equation (1) and its generaliza- 
tions with the use the Crocco type transformation 

Reduction of order of equation (1). Let us introduce the notation 

oz az^ 

Transferring the term mF'^ (here and below, the brief record of derivatives is used) 
to the right-hand side of equation (1), dividing both the resulting equation by F^- = <I>, 
differentiating with respect to z, and taking into account formulas (7), we eventually 
obtain the following equation: 

_ , ^ 9 v^z + rmf + q{t)r^ + pit) 

Now let us pass in equation (8) from the old variables t, x, F = F{x, t) to the new 
variables t,ri,^ ^ rj), where 7] and <i> are defined by formulas (7). The derivatives 
are transformed as follows: 

d drj d TP ^ ^ 

z oz or] OT] orj 

d d di] d d ^ d 
dt dt dt drj dt ^* drj 



As a result, equation (8) reduces to the second-order equation 

?/ did fmrf+qri+p\ 9^$ 

which can be rewritten as follows: 

— + [rmf +qrj+p)— = [(2m - 1)t] + q]^ + — . (9) 

Here and below, for brevity, the arguments of the functions p{t), q{t), and s{t) will 
often be omitted. 

Note that, in the degenerate case (inviscid fluid, — 0), the initial nonlinear 
second-order equation (1) is reduced to the first-order linear equation (9), which can be 
completely integrated using the method of characteristics. 

If a solution of the initial equation (1) is known, then formulas (7) define a solution 
to equation (9) in the parametric form. 

Let $ = $(77, t) be some solution to equation (9). Then, the corresponding solution 
to the initial equation (1) can also be presented in parametric form as follows: 

ds , . ^ f sds ,,, . 

where tp{t) is an arbitrary function (in the integrals, t is treated as a parameter). 
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Transformation of system (1), (2). Assuming that F^z ^ 0, let us pass in system 
(1), (2) from the old variables t, z, F to the new variables t, 77, $ according to formulas 
(7). Then, equation (1) is transformed into equation (9) and equation (2) is transformed 
into the following equation: 

dG d^G 
qT]+p)- k7]G = ly^^-^ + s(t). (10) 



dG , 2 



In deriving this equation, we used a representation for the mixed derivative obtained 
from equation (1). 

At m = fc, equation (10) has exact solutions of the following form: 

G = At? + S$ + C, (11) 

where A = A(t), B = B{t), and C ^ C{t) are unknown functions determined by an 
appropriate system of ordinary differential equations (at m = k ^ 1, we have B ~ 0). 
To prove this fact, one should substitute expression (11) into equation (10) and take 
into account equation (9). Formula (11) will be used below for the representation of 
solutions to equation (2) via solutions to equation (1). 

Some generalizations. Let us consider the nonlinear nth-order equation 



d'^F 



+ [a{t)F + b{t)z 



d^F 



H 



OF 
oz 



Q2p gip 



dz^ ' dz^ 



d"F 
9z" 



(12) 



dzdt dz^ 
which generalize equation (1) and allow for order reduction. Passing in equation (12) 
from t, X, F = F{x,t) to the new variables t, rj, ^ — ^{t,ri), where rj and $ are 
defined by formulas (7), we obtain the (n — Z)th-order equation 



a{t)r] + b{t) 



d_l_ 



d_ 



drj ' 



9"- 



d2 



(13) 



in which the high-order derivatives are calculated using the formulas 



d'^F d 



fc-2 



dz^ 



dz^ 



d_ 

dz 



d 



drj dz^ ^ ' dz drj^ 

In the special case of a second-order equation with n = 1, a{t) — —1, b{t) — 0, and 
H — H{Fz), equation (12) reduces to the Calogero equation, which was considered 
in [8, 12, 13]. It is evident that a more general equation (12) with n ^ 1, H = 
H{t, Fz), and arbitrary functions a{t) and b{t) (it is logical to call this the generalized 
Calogero equation), can be reduced using the Crocco type transformation to the first- 
order equation (13), which becomes linear after the substitution $ = l/^P. 



3 Representation of solutions to equation (2) via solu- 
tions to equation (1) 

The case of m = fc = 1. Let F — F{z,t) be a solution to equation (1). Then, by 
virtue of relations (11) and (7), equation (2) has the solution 

dF d^F 

G = A[+Aq + A— + B — , (14) 
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where the functions A = A{t) and B — Bit) satisfy the ordinary differential equations 

^;t + qA + {P + qt)A = s, (15) 
B[ + qB = 0. (16) 

The proof is carried out by eliminating the function G from (2) and (14), followed 
by the comparison of the obtained expression to both equation (1) and the equation 
resulting from the differentiation of equation (1) with respect to z. 

The general solution to equation ( 1 6) is as follows: B = C cxp J q dt^ , where 

C is an arbitrary constant. 

The case ofm — k^ 1. In this case, by virtue of relations (1 1) and (7), equation 
(2) has the solution 

1 dF 
G = -{A[ + Aq) + A—, (17) 
m oz 

where function A = A{t) satisfies the ordinary differential equation 

Att + qAf + {mp + qf)A:^ms. (18) 

The exact solutions to equation (1) can be found in [5, 8, 10, 11]. Formulas (14)- 
(18) are used to obtain the corresponding exact solutions to equation (2). 

4 Use of linear transformations for constructing exact 
solutions to equations (1) 

Using the following linear transformation with respect to the unknown function, 

F ^a{t)fiT,0+bit)z + c{t), 

r (19) 

^ =X{t)z + a{t), T= X'^(t)dt + Co, 



where a — a{t), b — b{t), c — c{t), A — X{t), and a — a{t) are arbitrary functions, 
we write equation (1) in the following form: 

^ + P(r)/ + 6(r)e + c(r)]^-mS(r)(^) = + (20) 

where 



q = —irlaqX + 2mabX — (aX)'f], p = -^{p + bq + mb^ — b'f). 



(21) 



The argument of the functions on the left-hand side of the equations is r and that on 
the right-hand side is t; the variables r and t are linked by the last relation in (19). 
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The presence of a large number (from five to seven) of arbitrary functions in equa- 
tions (1) and (19) allows us to construct various exact solutions to equation (1). 

Example 1. Assuming sequentially that (20) / — / = Ae^ + Be^^, f = 
^sin^ + B cosS,, and / — tanh^ at m = 1 in transformation (19) and equation (20), 
we can obtain the solutions given in [8, 11] by selecting appropriate arbitrary functions. 
Equation (20) is also satisfied by the functions / = (1 ± e^)^^ and / — tan(^ + A), 
which give new solutions. 

Example 2. Assuming that 

6 = C2, c^Cs, q^d, p^C^, (22) 

where C,i are arbitrary constants, we obtain an ordinary differential equation for the 
function / = /(^) from equation (20). In this case, relations (21) under conditions 
(22) represent a system of ordinary differential equations for the functional coefficients 
of transformation (19). In equation (1) with m = 1 and in transformation (19) with 
constraints (21) and (22), two functions can be set arbitrarily (it should be recalled that 
p = p{t) and q = q{t) are arbitrary functions). Stationary solutions / = /(^) of 
equation (20) generate nonstationary traveling-wave solutions (19) of the initial equa- 
tion (1). 

Example 3. Now we assume that 
fc+i 



ClT 2-, b = C2T~\ C^CzT 2, 



(23) 



where k and C„ are arbitrary constants. In this case, equation (20) admits self-similar 
solutions of the following form; 

where the function h — h{C,) satisfies the ordinary differential equation 

- C4) /4 + [C,h +(C2-^)C + Cs] h'l^ - mC\{h'^f ^ uh'l'^^ + C5. (24) 

Substituting expressions (23) into (21), we obtain a system of integro-differential 
equations for the functional parameters of the initial equation (1) and transformation 
(19). Note that, using the substitution A — y^^, and taking into account the relation 
T = ip{t) + Co, we arrive at a standard system of ordinary differential equations. As a 
result, a non-self-similar solution of the form (19) is obtained. 



5 Nonlinear analysis of stability/instability of solutions 

Analysis of stability/instability of solutions based on equation (2). Consider system 
(1), (2) with m — k — 1 and 5 = obtained in [11]. To analyze stability/instability 
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of solutions, we use formula (14) and equations (15) and (16) relating the solutions of 
system (1), (2). It is important to note that, in many cases, there is no necessity to know 
the explicit form of the function F. 

First, let us study problems with a stationary longitudinal velocity, which corre- 
spond to the case of F = F{z), p = const, and q = const. In this case, the solution to 
equation (15) depends on the sign of the discriminant A = 



A{t) = 



cxp 



cxp 



exp 




2 exp 



C2 cos 



\iClt + C2) 




(25) 



where Ci and C2 are arbitrary constants. Further, for the sake of simplicity, we assume 
i? = in equations (14) and (16). In the analysis, we consider the following two cases. 

1°. Nondegenerate case ^ 0. For g < (and arbitrary p) or p < (and 
arbitrary q), solutions (14) and (25) (with Ci 7^ 0) increase exponentially as t ^ 00. 
Therefore, the specified values of the parameters p and q determine the domain of non- 
linear instability of system (1), (2) for any limited stationary profile of the longitudinal 
component of the velocity F{z) (other than a constant). The point p ^ q = also 
belongs to the domain of instability of system (1), (2). 

Indeed, by choosing suitable values of the constants Ci and C2 and using equations 
(14) and (25), we can make the initial value of |G|t=o (interpreted as the initial pertur- 
bation relative to the trivial solution G = of equation (2)) smaller than any preset e. 
However, for q < (and arbitrary p) or p < (and arbitrary q), we have \G\ — *■ 00 as 
t — > 00. This means that arbitrary small perturbations of the solutions to system (1), 
(2) exhibit unbounded growth with time. 

Remark. If F ^ Fi at z ^ —00 and F ^ F2 as 2; ^ +00 (Fi, F2 = const), then 
solution (14) at A = 0, i? 7^ tends to zero as z ^ ±00. 

For q = and p > 0, solution (25) and hence solution (14) are periodic. The 
inequalities q > 0, p > {\p\ + \q\ 0) determine the domain of conditional stabiUty 
of the solutions under consideration. 

It is important to emphasize that (i) here we are speaking of the nonlinear insta- 
bility; (ii) all the results and solutions obtained above are exact (rather than linearized, 
which is the case in the theory of linear stability); and (iii) various assumptions, expan- 
sions, and approximations inherent in many nonlinear theories are not used either [2, 
14, 15]). 

Example 1. The stationary spatially periodic solution 

F = asin(crz + b), G = [p = —a^cP'^ q = va^^ 
of system (1), (2) is unstable for any values of a, 6, and cr (a 7^ 0, <t 7^ 0). 
Example 2. The stationary monotonic restricted solution 

F -6i^crtanh((Tz + 6), G = (p = 0, q ^ Si^a"^) 



of system (1), (2) is stable. 
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All conclusions concerning the stability/instability of solutions presented above, 
as well as formulas (14) and (25), remain valid for any nonstationary solutions F = 
F{z, t), G — G{z, t) (under the condition that the derivative F^ ^ is bounded) of 
system (1), (2) with p = const and q = const. 

According to the above considerations, three quarters of the plane of parameters p, 
q correspond to non-stationary solutions. It is important to note that the flow instability 
described above is not associated with a specific velocity profile and is realized due to 
equation (2) governing the transverse components of the fluid velocity. Since the fluid 
viscosity v does not enter equations (15) and formulas (25), the above results do not 
depend on the Reynolds number; i.e., the instability of solutions takes place not only at 
large but also at small Reynolds numbers (0 < Rc < cxi). 

Remark. Likewise, we can use equations (14) and (15) to study instability of the 
nonstationary solutions of system (1), (2) with variables p = p{t) and q = q{t). 

2°. Degenerate case F^ = 0. Let 

F = a = const (if p = 0). (26) 

Then any solution to equation (2) that, when passing from t, z to the new variables 
^ = z — at, reduces to the classical heat equation, is stable for any values of the 
parameters a and q. 

Analysis of stability of solution F — const of equation (1) for p = 0. 

1°. Let us study stability of the trivial solution = of equation (1) forp = and 
various values of the parameter q — const. Equation (1) admits the following exact 
solution; 

F = £6**^^+^*, \ = q- vk^, (27) 

where e, k, and A are real quantities. This solution is also a solution to the linearized 
equation (1) with the quadratic terms discarded. The absolute value of the difference 
between solution (27) and the trivial solution at the initial moment of time is |e| (this 
difference can be made arbitrarily small by choosing a suitable e value). 

For q — vk^ > 0, the trivial solution will be unstable, whereas it will be stable for 
q — vk^ < 0. The stability boundary is the parabola q — vk"^ on the k, q plane. For 
decreasing fluid viscosity, v ^ Q (which corresponds to increasing Reynolds numbers), 
the branches of this parabola tend to the line g = and the domain of instability 
expands to become, in the limiting case, the entire upper half-plane q > 0. Increasing 
V ox k results in the expansion of the stability domain. Since the parameter k can be 
set arbitrarily, then, for any g > we can achieve instability of the trivial solution by 
choosing an appropriate k value. 

2°. Consider an arbitrary stationary solution of the form (26). Instead of solution 
(27), we use the following function: 

F = ^e^H^-at)+\t X = q- i/fc2, (28) 

Owing to property (3) with — —at, this function is also a solution to equation (1). 
The absolute value of the difference between equations (26) and (28) at the initial time 
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can be made arbitrarily small by choosing a suitable e value. All criteria of stability 
and instability of solution (26) depending on the parameters k and q remain the same 
as those for the trivial solution. 

Remark. It follows from the above results that, forp = and q < 0, only a constant 
profile of the longitudinal velocity component {F = const) is stable. 



6 Transformation and exact solutions of the boundary 
layer equations 

Nonstationary equations of the plane boundary layer in terms of the stream function w 
are reduced to a single third-order equation [8]: 

d^w dw d^w dw d'^w d^w 

dtdy dy dxdy dx dy'^ dy^ ' 

Let us pass from t, x, y, and w to the new variables t, x, rj, and $, where ?y, $ are 
the generalized Crocco variables defined as follows: 

dw ^, . d'^w 
r, = —, $(,,i,^) = _. 

As a result, we obtain the following second-order equation, 

which is reduced by the substitution <i> = 1 to the nonlinear heat equation: 

9f + 9^ + ~ ""a^ ■ 

1°. Let us first consider the special case of f{x, t) — f{t) and for special solutions 
of the following form: 

«' = Z(C,r), ^^x-Tjt + jtf{t)dt, r = ^t^ 

Here, we have the solvable equation 

which can be reduced to the linear heat equation [8]. 

2°. Consider a more general case of f{x, t) ~ f{t)x + g{t) and seek the following 
special solutions: 

* = Z{t r), C = v{t)x + ^(t)77 + e{t), r = / ^^^(i) dt, 
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where the functions = ip{t), ip = and 9 = 9{t) are determined by the following 
linear system of ordinary differential equations 

(^; + /V;-o, ^^ + ^ = 0, 0; + .gv = o. 

As a result, we also arrive at the solvable equation (29). 
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